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ABSTRACT 

Modern data of the extinction curve from ultraviolet to near infrared are revisited to study the 
property of dust grains in the Milky Way (MW) and the Small Magellanic Cloud (SMC) . We confirmed 
that the graphite-silicate mixture of grains yields the observed extinction curve with the simple power- 
law distribution of the grain size but with a cutoff at some maximal size: the parameters are tightly 
constrained to be g = 3.5 ± 0.2 for the size distribution and the maximum radius flmax = 0.24 ± 
0.05 Ilia, for both MW and SMC. The abundance of grains, and hence the elemental abundance, is 
constrained from the reddening versus hydrogen column density, E{B — y)/A^H- If we take the solar 
elemental abundance as a standard for the MW, >56 % of carbon should be in graphite dust, while 
it is <40 % in the SMC. This disparity and the relative abundance explain the difference of the two 
curves. We find that 50-60 % of carbon may not necessarily be in graphite but in the amorphous or 
the glassy phase, or up to 15 % in silicon carbide. Iron may also be in the metallic phase or up to ~80 
% in magnetite rather than in silicates, so that the Mg/Fe ratio in astronomical olivine is arbitrary. 
With these substitutions the parameters of the grain size remain unchanged. The mass density of 
dust grains relative to hydrogen is pdust/pH = 1/(120+}°) for the MW and 1/(760+™) for the SMC 
under some abundance constraints. We underline the importance of the wavelength-dependence slope 
of the extinction curve in the near infrared in constructing the dust model: if oc A"'' with 7 ~ 1.6, 
the power-law grain-size model fails, whereas it works if 7 ~ 1.8-2.0. 

Subject headings: dust, extinction - ISM: abundances - ISM: general - Galaxy: abundances 



1. INTRODUCTION 

Submicron-size grains cause attenuation of light from 
ultraviolet (UV) to near-infrared (NIR) wavelengths. 

' Larger grains may form small astronomical objects. 
Thirty years ago, Mathis, Rumpl, & Nordsieck (1977, 
hereinafter MRN) showed that dust composed of a mix- 

I ture of silicate and carbonaceous grains accounts for the 
extinction curve from the UV to optical wavelengths, 
from which they derived the size distribution of grains 

I for a = 0.005 to 0.25 fim. Only with 9 bins of the his- 
togram, they advocated that the size distribution is con- 
sistent with the power-law a~'^ with the index q from 3.3 

[ to 3.6."* This dust model is consistent with the fact that 
Mg, Si, and Fe are highly depleted in diffuse interstellar 
matter. These three elements arc also major constituents 

I of astronomical silicate with the corroborating fact that 
the abundance of dust needed to account for the observed 
extinction is on the order of magnitude of their cosmic 
abundance, and hence the three elements could domi- 
nantly be locked in dust. 

Independently, the argument appeared that small as- 
tronomical bodies in frequent collisions would obey a 
power-law distribution in size (Dohnanyi 1969; Hellyer 
1970; Biermann & Harwit 1980; Dorschner 1982). It was 
argued that the power-law derived by collisional equilib- 
rium has typically q ~ 3.5 (Dohnanyi 1969; Pan & Sari 
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power-law: it results in q = 3.8 ±0.7 for graphite and q = 3.5 ± 0.6 
for silicate, where the error stands for the dispersion. 



2005), which agrees with the power favoured by interstel- 
lar dust from extinction studies. For larger sizes the dis- 
tribution of small astronomical objects (> 1-10 km) does 
not conflict with this power, while the available piece of 
information is limited to derive it more accurately. 

These arguments indicate that the size distribution of 
grains may be a key to understanding the formation of 
grains and small astronomical objects. The size distri- 
bution of grains has been enquired in the literature from 
time to time (e.g., Draine & Lee 1984, hereafter DL84; 
Kim, Martin, & Hendry 1994; Weingartner & Draine 
2001, hereafter WDOl; Clayton et al. 2003a; Zubko, 
Dwek, & Arendt 2004) . Later studies are weighted more 
in favour of accurately reproducing the extinction curve 
that is progressively more accurately measured and is ex- 
tended to the NIR region, to find the elaborate size dis- 
tribution of grains. A significant variation has become 
apparent in the extinction curve depending on the line 
of sight (Mathis & CardeUi 1992; Clayton et al. 2000). 
The variation typically amounts to Ry = 2.2-5.5, if the 
one-parameter formula of CardeUi, Clayton, & Mathis 
(1989, hereafter CCM: their formula is referred to as the 
CCM curve) is adopted, who showed that the family of 
the extinction curve is well described by the formula with 
the single free parameter, Ry = Ay / E{B — V), where 
Ay is the extinction in the V band and E{B — V) is the 
reddening. It should be asked how the variation of the 
extinction curve is translated to that of the property of 
dust grains. It is also a problem to be asked what can be 
constituents of grains other than silicate and graphite. 

Further to the shape of the extinction curve, the size of 
the extinction, usually represented with E{B — V) is the 
quantity that gives an important constraint on the abun- 
dance of dust grains. An important question is whether 
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the grain parameters derived from the extinetion curve 
are consistent with those from the size of extinction, and 
whether they are consistent with the elemental abun- 
dance. 

It has been known that the extinction curve for the 
Small Magellanic Cloud (SMC) is markedly different 
from that for the Milky Way (MW) in that it lacks the 
feature at 2175 A and it shows a significantly steeper rise 
to the far UV side beyond 2000 A. What would cause this 
difference is the problem to be asked. The SMC type ex- 
tinction curve is also indicated for interstellar clouds such 
as Mg II absorbers (York et al. 2006); it is shown that the 
majority, typically '^70 %, of the cloud obeys the SMC 
type extinction law. 

In this paper we revisit studying gross but generic fea- 
tures of dust grains using the modern data of extinction 
curves from 1150 A to 2.2 /xm derived for the MW and, 
in the other extreme, the SMC, along with modern op- 
tical data of the relevant dielectric material. We exam- 
ine whether some departure from power-law size distri- 
butions, other than cutoffs in the size, is compelling in 
reproducing the observed extinction curve. We are also 
interested in how dust in the SMC should differ from that 
in the MW in its properties. Furthermore, we would like 
to see whether the graphite-silicate model is unique and 
what is the possible range of dust to hydrogen mass ratio. 
We assume spherical grains and their size distribution 
obeying a power-law allowing for a truncation at some 
maximum size. The power-law distribution, at least, can 
be understood with simple physics for the evolution of 
grains. We do not treat polarisation and infrared emis- 
sion from dust since either requires the knowledge other 
than the property of dust grains and so, in turn, requires 
extra assumptions and introduces more uncertainties. 

In Section 2, we review the data of extinction curves 
used in the present study, and we note a problem in the 
average NIR extinction curve. After defining the model 
of interstellar dust used in this paper in Section 3, we 
search, in Section 4, for the dust model that could repro- 
duce the observation, present parameters of the model 
and ask various possibilities of the favourable compo- 
sition of grain species from the extinction data for the 
MW. We carry out a similar analysis for SMC dust in 
Section 5. Section 6 gives a summary of this study with 
discussion. 

2. EXTINCTION CURVES 

2.1. Milky Way extinction 

We take the extinction curves derived from UV to NIR 
for 328 stars in Fitzpatrick & Massa (2007, hereafter 
FM07) with the aid of stellar atmosphere, replacing the 
traditional method using rcddened-comparison pairs of 
stars. The extinction curve ranges from 1150 A to 2.2 
/xm. We arc primarily interested in the global average 
of the extinction curves, so we consider a bunch of the 
curves at several specific reference wavelengths, at which 
we derive the allowed ranges of extinction. We take 9 ref- 
erence wavelengths in total with the V band to give the 
normalisation. We consider wavelengths corresponding 
to UBJHK and three wavelengths, A = 0.2175 fim, 0.16 
/im, and 0.125 /xm in UV, which characterise the hump, 
the bottom in the UV region, and an arbitrary chosen 
wavelength in the rise of the Galactic extinction curve 



TABLE 1 

Interstellar Extinction Data for the Milky 
Way 



Reference Wavelength 


1/A 




v)obs 






(fim 1) 


lower 


upper 


0.125 


(far-UV rise) 


8.00 


2.573 


3.894 


0.16 


(far-UV dip) 


6.25 


2.167 


3.003 


0.2175 


(UV bump) 


4.60 


2.712 


3.625 


0.36 


( U band) 


2.78 


1.449 


1.62 


0.44 


(B band) 


2.27 


1.252 


1.331 


1.26 


(J band) 


0.80 


0.211 


0.278 


1.62 


(H band) 


0.61 


0.127 


0.166 


2.15 


{K band) 


0.46 


0.077 


0.101 


Ran 


ges from the CCM formul 


a for NIR 


1.26 


(J band) 


0.80 


0.267 


0.299 


1.62 


(H band) 


0.61 


0.171 


0.191 


2.15 


(K band) 


0.46 


0.110 


0.123 



Note. — The allowed ranges of the interstellar 
extinction {Ax/Ay)^i^i^ at the reference wavelengths 
constructed from the data of extinction curves of 
FM07. The ranges of the NIR extinction from the 
CCM formula is for Ry = 2.75-3.60. In optical re- 
gions, the ranges virtually agree with the 1 a ranges 
of FM07. 

towards shorter wavelengths, respectively. We take the 
l-cr allowed ranges corresponding to 68 % (224 curves) of 
the 328 curves at each reference point. Our wavelength 
mesh is too coarse to study the 2175 A feature; so it is 
considered separately. We do not treat diffuse interstel- 
lar bands which are not apparent in the FM07 extinction 
curve. Possible line features beyond the K band arc not 
treated. 

The l-cr ranges of A\/Av are shown in Table 1 and 
in Figure 1: (a) for the entire wavelength range consid- 
ered and (b) for NIR in an expanded scale. Figure 1(a) 
shows that the CCM curve with Ry = 2.75-3.60 and the 
dust-model calculation of WDOl (which gives Ry = 3.1) 
are both consistent with the data of FM07 in the optical 
and the UV regions. In the NIR region, however, the ex- 
panded figure. Figure 1(b), indicates that both curves arc 
off from the FM07 data by Ril a. This arises from the fact 
that CCM considered the NIR extinction that follows the 
power-law A\/Av oc with the index 7 = 1.61 (Rieke 
& Lebofsky 1985), whereas the FM07 data are consistent 
with a steeper index, 7 — 1.84. The steeper index (Mar- 
tin & Whittet 1990) has gained more supports in the 
recent work: for example, Fitzpatrick & Massa (2009) 
reported the index 1.78-2.0 and that it varies according 
to the line of sight. 

We take the FM07 data as our prime choice, while 
tentatively retaining the possibility that the NIR power 
index is moderate, as with CCM and WDOl. We also 
consider this possibility and study the implication on the 
dust model. We take the range of the CCM curves with 
Rv ~ 2.75-3.60, which correspond to 1 ct of the extinc- 
tion data of FM07 in the UV and optical range, as seen 
in Figure 1(a). Our 1 sigma range for NIR covers the 
variation of the index reported by Fitzpatrick & Massa 
(2009), while CCM and WDOl are beyond 1 a. The 
difference in NIR leads to a significant difference in the 
conclusion concerning the model. 

In addition, we consider the absolute size of extinction. 
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Fig. 1. — (a) Interstellar extinction curves from the 328 stars in 
the Milky Way from FM07. Our 1 cr (thick solid bars) estimates 
are denoted at our reference wavelengths. The plus symbol shows 
the V band used as the normalisation, (b) The expanded figure 
for the extinction in NIR. The bars are 1 a. The extinction curves 
of the CCM formula with Ry = 2.75 [the upper dashed curve in 
(a) and the lower dashed curve in (b)] and Ry = 3.60 [the lower 
dashed curve in (a) and the upper dashed in (b)], and the model 
of WDOl {Rv = 3.1, the solid curves) arc shown. 

or reddening. Bohlin, Savage, & Drake (1978) obtained 
using 75 data points Nii/E{B - V) = 5.8 x 10^^ cm-^ 
mag~^, which is widely adopted in the literature, and 
claimed that the data for different lines of sight rarely 
fall beyond the lines 1.5 or 1/1.5 times the value indi- 
cated in this expression. Actually, we see in their figure 
(their Figure 2b) that about 85 % of stars are located 
between the two lines. (Most of the deviants are towards 
the smaller TVh side.) The Nh value includes a 25 % con- 
tribution from H2 molecules. Unfortunately, they have 
not given the error or the dispersion of the fit. So we have 
re-fitted their data to obtain the 1 a error. We have tried 
to include the recent enlarged data set of Gudennavar et 
al. (2012). After the selection similar to that in Bohlin 
et al. (1978) we obtain, using 174 data. 



Nil/ E{B - V)^ {5.7 ± 1.7) X 10^^ cm- 



mag \(1) 



where the error stands for the dispersion of the fit. This 
is in good agreement with the original Bohlin et al. 
(1978) result, including the size of the dispersion we re- 
estimated from their data (±1.7). 



TABLE 2 

Extinction Data in the Small Magellanic 
Cloud 

Reference Wavelength 1/A {A\/-Av)oha 
(fim) {fim~^) upper lower 



0.125 


(far-UV rise) 


8.00 


6.752 


4.839 


0.16 


(far-UV dip) 


6.25 


5.005 


3.813 


0.2175 


(UV Bump) 


4.60 


3.946 


2.772 


0.36 


( U band) 


2.78 


1.892 


1.521 


0.44 


(B band) 


2.27 


1.488 


1.303 


1.26 


(J band) 


0.80 


0.324 


0.108 


1.62 


(H band) 


0.61 


0.184 


0.000 


2.15 


{K band) 


0.46 


0.060 


0.000 



2.2. Extinction in the Small Magellanic Cloud 

Only a handful of sightlines are studied for the extinc- 
tion curve for the SMC. Gordon et al. (2003, hereafter 
G03) presented the curve towards 5 stars, AzV 18, AzV 
23, AzV 214, AzV 398, and AzV 456 using the traditional 
reddened-comparison pair of stars. More recent result of 
Cartledge et al. (2005) is consistent with G03. 

As known, the bump at 2175 A is not apparent except 
for the sightline towards AzV 456, which is located in a 
'wing' region of quiescent star formation, unlike the other 
four that pass through the star- forming bar of the SMC. 
We include the extinction curve towards AzV 456 in our 
consideration. 

For our measure we define, because of the paucity of 
data, the '1 a ranges' of the extinction to be the maxi- 
mum and the minimum of the five curves at each refer- 
ence wavelength. If the lower value becomes negative, it 
is set to zero. The ranges of extinction are given in Table 
2 and plotted in Figure 2. 

The recent analysis for the column density of neutral 
hydrogen A'hi is by Welty et al. (2012). The abundance 
of H2 molecules is estimated to be 2 % of hydrogen, which 
is compared to 25 % for the MW (Bohhn et al. 1978), 
giving an example of molecular formation depending on 
environment. Their estimate is consistent with that of 



Tumlinson et al. (2002)^ 27Vh2/A^h 
et al. (2012) have given 



1 



-1-0.5 
-0.2 



%. Welty 



N^/E{B-V) = 2.3_^^ X 10 



99 —9 —1 

cm ^ mag . 



(2) 



This is larger than the MW value given in Equation (1) 
by a factor of 4.0 ± 2.4. Russell & Dopita (1992) give 
a summary of the elemental abundance for the SMC, 
which indicates that its metallicity is 5.6 times smaller 
than the solar. The inferred dust abundance of the SMC, 
relative to the MW, is consistent with, or somewhat less 
suppressed compared with, the inferred heavy element 
abundance. 

3. DUST MODEL 

With spherical particles uniformly distributed in inter- 
stellar space the total extinction, A\, along the line of 

^ The estimate of Cartledge et al. (2005) for H2 for some of the 
5 stars wildly varies from 5 % to 50 %. We quote that the cosmic 
global abundance of 2Af(H2)/Af(H) is inferred to be 0.30 (Fukugita 
2011). 
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Fig. 2. — (a) The extinction curves for the SMC derived from 
five stars by G03, with 1 a (thick soUd bars) significance denoted 
at the reference wavelengths. The plus symbol corresponds to the 
V band, (b) The expanded figure for the NIR extinction. The 
extinction curve given by WDOl (for AzV 398) is shown by the solid 
line. The extinction curve for the MW, taken from Figure 3 below 
(our model), is also added by thin dotted curves for comparison. 

sight I, is 



A, 



1.086 



dl 



(3) 



where (5^'^*(a) is the extinction efficiency defined as the 
ratio of the extinction cross section (T\j{a) at wave- 
lengths A to the geometric cross section na^ for grain 
species j and is calculated using Mie scattering with the 
laboratory optical data of dielectric constants: nj{a)da 
is the number density of the grain of species j with radii 
between a and a -\- da. 

The "standard composition" of dust is graphite and 
silicate. Following DL84 wc take "astronomical silicate" 
(we simply call it silicate unless otherwise stated) with 
the composition MgFeSi04. We assume the simple ap- 
proximation that dust grains are bare refractory particles 
without substructure. For graphite, we calculate ^A^gra 
using the usually adopted — § approximation for the di- 
electric constant components perpendicular and parallel 
to the basal plane (e.g., DL84) to represent anisotropy. 

In order to see their possible importance we also con- 
sider other refractory components and carbonaceous ma- 



terials, which arc likely to condense into grains, in so far 
as optical properties are known. We consider 10 grain 
species, graphite (gra; optical data from Draine 2003b), 
glassy carbon (glC; Edoh 1983), amorphous carbon 
(amC; Zubko et al. 1996), silicon carbide (SiC; Choyke 
& Palik 1985), astronomical silicate (asil; Draine 2003b) 
of the chemical composition of olivine (Mga:Fei_aSi04) 
with a; w 1, forsterite (Mg2Si04; Semenov et al. 2003), 
pure iron (Fe, Semenov et al. 2003), magnetite (Fe304; 
Triaud)^, troilitc (FcS; Semenov et al. 2003), and corun- 
dum (AI2O3; Toon et al. 1976). We do not take enstatite 
(MgSiOa). One may expect that the optical property of 
MgSiOa is not much different from that of olivine, and its 
contribution may not disturb much the calculation with 
astronomical silicate.^ We ignore Si02, since its Q factor 
is small. Inclusion of Si02 only increases the amount of 
Si locked in dust grains. The Ca- and Ti-bearing grains 
may constitute some components of dust, but they are 
not considered because of their small cosmic abundance. 

We do not include polycyclic aromatic hydrocarbons 
(PAHs) separately. The bump at 2175 A in the Galac- 
tic extinction curve can be accounted for either with or 
without PAH in so far as small graphite grains (molecules 
with the effective radius of <200 A) are included (Stecher 
& Bonn 1965; Joblin et al. 1992; Clayton et al. 2003b). 

Dust grains dominated by a single size do not give the 
observationally obtained extinction curve. Some distri- 
bution over the size is necessary. A typical distribution 
that is known to work is the power law, which we take 
here: 



nj{a) = nHi^j(a/ao) 



-9j 



(4) 



where nn is the hydrogen number density, Kj is the frac- 
tion of species j, and ag is a size of the normalisation. We 



limit the range to 



< a < a-n 



We assume the 



same grain-size distribution independent of grain species. 



We take qj, 



•''maxj 1 



and the condensed fraction /i_ 



3 1 



i.e., the fraction of element i contained in grain species 
J, as parameters, and find a set of the parameters that 
satisfy l-cr of extinction data at all our reference wave- 
lengths. We have studied the dependence on the param- 
eter a„iinj , but we have seen that it is small in so far 
as we take it to be smaller than 0.005 /im. The min- 
imum size flmin affects the extinction curve in the way 
that increasing it diminishes the rise in the far UV for 
wavelengths roughly shorter than 0(27r) x Omin- With 
our choice of amin = 50 A it affects little the curve for 
A > 1150-1250 A we consider. The increasing variation 
of amin also affects the 2175 A bump caused by graphite, 
reducing the hump, but the resulting change of the ex- 
tinction curve from the decreasing variation of Omin from 
0.005 /xm to zero is small, typically <l/3 the 1 a error 
range. This parameter is not well determined if it is left 
as a free parameter, unless the extinction data of shorter 
wavelengths are used. So we fix amin at 0.005 /im. The 
change of lower cutoff affects little the determination of 
other parameters, unless g > 4 that is deep in the region 



.astro.uni-jena.de/Laboratory/OCDB/ 1 



A. Triaud, [http://v 
oxsul.html 

^ The imaginary part of the presently available refractive index 
of enstatite (Dorschner et al. 1995) deviates largely from that of 
forsterite in the UV region, while the real part differs little. The 
Q factor of enstatite differs little from that of forsterite. So we do 
not treat enstatite separately in this paper. 
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not allowed in our study. It increases the total amount 
of mass, say, by 10 % if Cmin is reduced from 0.005 fim 
to 0.001 /im. A possible increase of the lower cutoff for 
graphite is considered for SMC dust separately, however. 
The extinction curve we refer to is Ax/ Ay normalised in 
the V band: 



Av 



E , K,A 



(5) 



where j is the component contribution to Equation 
(3). 

We also calculate 



E{B - V) 



= 1.086 ^A'j 



da{aB,j 



'^vj) — 



-i3 



(a) 



- q=3.5 
=0.24 /um 



C.gra/^Fe.asil" 



=0.7 



1=3.2 
1=3.8 

0.34 ixm 
0.2 



max 
fc,gra/f, 



Fe.asil 



(6) 



The size of reddening per hydrogen leads to the abun- 
dance of dust that can be compared directly with esti- 
mates of the elemental abundance. Assuming that the 
abundance does not vary from place to place, we take 
as a standard for the MW the elemental abundance es- 
timated by Grevesse & Sauval (1998) (hereafter 0898).^ 
Asplund et al. (2009), using the 3D calculation, claimed 
the solar abundance generally lower than that of GS98 
by as much as 30 % especially for C, O and a few others. 
A comparable 3D calculation (Caffau et al. 2011), how- 
ever, has given the abundance that is lower than GS98 
only by 10 %. 

The elemental abundance is generally not tightly con- 
verged among the authors. For instance, for the car- 
bon abundance, which is of one of our major concerns, 
the GS98 value C/H (3.3 ± 0.5) x 10"'^ may be com- 
pared with 3.62 X lO"'' of Anders & Grevesse (1989), 
2.69 X 10-4 of Asplund et al. (2009), (1.9-2.9) x lO"'' of 
Cardelli et al. (1996), 2.14 x 10""* of Nicva & PrzybiUa 
(2012), 2.45 x 10""* of Lodders (2010), and so forth. The 
iron abundance, for which GS98 gives log(Fe/II) -|~ 12 = 
7.50 ± 0.05, too varies between 7.45 and 7.66. Keeping 
these uncertainties, in mind, we take as our default GS98 
that leads to the structure of the sun in a satisfactory 
manner. 

For the SMC we take the composition given by Russell 
& Dopita (1992), with which total metallicity is 1/5.6 
times the solar. The abundance of refractory elements, 
Mg, Si, and Fe, is smaller than their MW values by sup- 
pression factors, 1/3.3 to 1/4.6. The abundance of O 
and C in the SMC seems relatively more strongly de- 
pressed, by 1/(6.2-6.3) times, compared with the solar 
abundances. This leads to the significant result that the 
ratio of carbonaceous material to silicate in the SMC is 
50 % lower than the corresponding value for the MW. 

4. RESULTS FOR MW DUST 

4.1. Single grain species 

We first study the extinction curve with a single grain 
species. We consider the allowed region in q and Omax 
plane, with which the dust model satisfies the 1 a ranges 
of FM07 extinction curve. We find that there are no 
overlaps among the three regions derived from the UV 
group (0.125 /xm, 0.16 /xm, and 0.2175 /im), the UBV 

* We remark that GS98 agrees with the earUor table of Anders 
&; Grevesse (1989) up to oxygen, for which GS98 value is lower by 
0.1 dex. 
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Fig. 3. — (a) Typical extinction curves from the graphite-silicate 
models, and (b) the expanded figure of (a) for NIR wavelengths. 
The error bars stand for the observed 1 a ranges. The thick solid 
curve is with q = 3.5, Omax = 0.24 /xm, and /c,gra//Fc,asil = 0.7, 
taken as our fiducial. The other curves denoted by thin curves are 
examples where one of the three parameters is shifted to demon- 
strate the response to the parameter. The legend shows the pa- 
rameters when changed from the fiducial choice. 

group, and the NIR group (J, H, and K) for any species 
of grains we considered. The regions required for UV 
and UBV groups are always disjoint. It often happens 
that no consistent parameters exist for J, i/, and K. We 
cannot make a dust model, which explains the extinction 
curve over the wide range of wavelength only with a sin- 
gle grain species. For instance, astronomical silicate gives 
to steep rise in the far UV if Cmax is chosen to account 
for the optical extinction curve, and the predicted NIR 
extinction is too small by a factor of 3 to 4. We see that 
these problems are readily offset by carbonaceous grains. 
This leads us to suppose that the mixture of these two 
species would give the correct extinction curve. 

4.2. Graphite- silicate dust model 

The two-component dust model, consisting of graphite 
and astronomical silicate, has been widely taken for in- 
terstellar dust since MRN. Mg, Si, and Fe are similar 
in the cosmic number abundance (within 20 %). If we 
assume that Si is all condensed into astronomical sili- 
cate, /Si.asil = 1, then /Mg,asil = 0.93 and /Fe.asil = 1, if 

Fe:Mg=l:l (see Drainc 2003a) with the GS98 abundance. 
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Here iron is slightly (10 %) deficient and extra Si and Mg 
may condense into forsterite. Oxygen locked in silicate 
is /o,asii = 0.26. Carbon is also depleted in interstellar 
matter. Following Sofia et al. (2011), it is approximately 
60-70 %, i.e., 

/Cgra = 0.6 - 0.7 (7) 

if all condensed carbon is in graphite. This is higher than 
their earlier estimate, 30-40 The depletion given by 
Cardelli et al. (1996) is also consistent with eq.Q if the 
GS98 abundance is adopted. 

The model attains good fits from far UV to NIR within 
1 (7 of the data with 

g = 3.5 ±0.2 and a,nax = 0.24+JJ-J,^. (8) 

—0.05 

The range of /c,gra//Fc,asii is 0.25-2.23, allowing for 
/Fe,asii < 1- Figure 3 presents examples of the result- 
ing extinction curve from our model at the central val- 
ues of eq.® with /c,gra//Fe,asii = 0.7. Wc draw sev- 
eral curves in addition that are somewhat away from the 
best fit value (say, about 1.5 a of the observation) to 
indicate how is the curve shifted with the variation of 
g, Omax, and /c,gra//FG,asii, the Other parameter that af- 
fects the shape of the extinction curve. The mass den- 
sity ratio of graphite to astronomical silicate is given by 
Pgra/Pasii = 0.73/c,gra//Fo,asii- Thc rolc of graphite is, 
in addition to producing 2175 A bump, to increase the 
extinction in NIR, which is too small with silicate alone 
when normalised in thc optical region. Graphite makes 
thc rise in the far UV more moderate. 

In Figure 4, we show in g — a„iax plane the region where 
the model gives the 1 a ranges of the extinction for the 
choices of /c,gra//Fe,asii = 0.7 and 0.2. The overlapping 
region is seen for UV, UBV, and NIR for /c,gra//Fo,asii = 
0.7, validating the graphite-silicate model. Such allowed 
regions disappear for /c,gra//Fo,asii < 0.25. 

Figure 5 shows the abundances of carbon and silicon 
in grains relative to hydrogen necessary to account for 
the extinction curve. The allowed region is in a belt 
/c,gra//Fc.asii =const, running from bottom left to top 
right. We also indicate the region allowed from E{B — 
V)/Nb_- This leads to the elemental abundance (carried 
by dust), 

log (C/H)-M2 = 8.4±0.3 and log (Si/H) + 12 = 7.6±0.4. 

(9) 

The shade shows various abundance estimates and thc 
dashed lines indicate GS98, 8.52 for C and 7.55 for 
Si. The ranges required from the interstellar extinction 
are consistent with other abundance estimates, includ- 
ing GS98, at least for C, Si, Mg, and Fe, we consider 
here. If we take the GS98 value for the total abundance 
of C and Si including the gas phase, the figure, showing 
log (C/H) + 12 > 8.27, means 

/Cgra > 0.56, (10) 

assuming that /si, grain = 1- This carbon fraction in 
graphite is consistent with the depletion estimated for 
interstellar matter referred to in eq.([7|)^°. 

^ Sofia ct al. (2004) used only weak absorption lines of car- 
bon, while their updated analysis (2011) includes strong absorption 
lines. 

The carbon abundance required from the extinction is still 
consistent with the lower abundance of Asplund et al. (2009), but 
then the depletion must be as high as >80 %. 
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Fig. 4. — Allowed parameter regions for graphite-silicate models 
in q and amax plane for the 1 cr range in thc UV group (blue), the 
UBV group (green), and the NIR (JHK) group (red). We take 
/c,gra//Fc asil = 0-7 in (a) and 0.2 in (b), the latter of which is 
a model out of 1 cr. Thin dotted curves show the parameters for 
thc 1 (T range in the NIR (JHK) group when thc NIR wavelength 
dependent power is 7 = 1.6 as with CCM or WDOl. 

A similar consideration is made for thc 1 a ranges 
from the CCM formula, where the NIR extinction takes 
a smaller power index of the wavelength dependence. 
The range of q and Omax required for the NIR extinc- 
tion curve is shown in Figure 4 above with dotted curves 
(those for UV and UBV remain unchanged). We see that 
the region allowed simultaneously for the three groups of 
color bands does not exist. The region for NIR is al- 
ways disjoint from those for UBV and UV. No graphite- 
silicate model is consistent with the observed extinction, 
if the NIR extinction power is ~A~^'®, for grains with 
the power-law size distribution. 

Thc steeper wavelength dependence for the NIR power 
thus looks more easily accommodated from the model 
point of view. We remark that WDOl tweaked signif- 
icantly the grain size distribution from thc power law, 
separately for graphite and silicate, adding different com- 
ponents so that the graphite-silicate model becomes con- 
sistent with the CCM-like extinction. This is also true for 
the size distributions obtained by other authors (listed 
earlier) when the CCM-like extinction is reproduced with 
the graphite-silicate model. Let us comment that thc so- 
lution of WDOl takes the Si, Mg, and Fe cosmic abun- 
dance larger by 30-50 % than the GS98 cosmic value. 
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Fig. 5. — The abundance of C and Si in dust grains relative 
to hydrogen for the MW to satisfy 1 sigma range of the extinction 
curve. The range within the oblique belt is allowed from the ex- 
tinction curve, and further restricted range in the middle is from 
E{B — V)/Nii. The two horizontal and vertical dashed lines show 
the total abundance of Si and C from GS98, with the shaded re- 
gions corresponding to the range of various estimates. 

We note that the aUowed regions in Figure 4 above 
stand for the ranges consistent with the variation of dust 
properties along hues of sight in the MW. They are well 
converged to narrow regions of q and a,„ax indicated in 
eq.®, despite an apparently significant variation of ex- 
tinction curves. The size property of grains is similar; 
only a small variation of the size parameter could cause 
the difference in the extinction curve. 

We summarise in Figure 6 the region allowed for Omax 
and q, marginalising over the ratio of condensed frac- 
tions of graphite to silicate for /c,gra//Fe,asii > 0.25. The 
curves are inlaid for /c.gra//Fc,asii = 0.4, 0.7, and 1. The 
1 a allowed range is realised with 3.2 < q < 3.7 and 0.19 
/im < Cmax < 0.34 /im. With increasing /c,gra, Omax be- 
comes smaller. The other set of curves is the constraint 
from E{B - V)/Nii for /c,gra = 0.4, 0.7, and 1 with 
/pe.asii = 1- Overlaps are seen between the two curves 
for /c,gra//Fo,asii = 0.7 and 1, but not for 0.4. 

For our model that satisfies 1 a constraints of both 
extinction curve and the E{B — V) size we have the 
mass density of dust relative to that of the hydrogen 
gas, Pdust/PH ~ l/(122]|^7g). The allowed range would 
be tightly constrained if the abundance is constrained to 
GS98 and further /Fe.asii = 1 is assumed: 

Pdust/PH = 1/(122 (11) 

The mass extinction constant i^oxt.A = ^A/^dust for the 
V band is 

K^^uv = (3.7 ± 0.5) X 10^ mag cm^ g"\ (12) 

which is compared to 2.8x 10'^ mag cm^ g~^ of WDOl. (If 
the abundance constraints are removed, eq. p^ becomes 
(3.6 ± 1.0) X 10^ mag cm^ g~^-) The ratio of the mass 
density of graphite to silicate is Pgra/pasii = 0.51 ^[J'^p, 
or in terms of the number of C and Fe atoms in grains, 
Nc/Nfo = 7.3+^4, corresponding to eg. pTj) . 



Interstellar Dust 7 




MW : 



6 



q 

Fig. 6. — Allowed regions of q and amax with which the model 
satisfies the 1 a ranges of the observed extinction for the MW, 
marginalised over the graphite-to-silicate ratios /c,gra//Fc,asil ^ 
0.25. The three contours inlaid in the allowed region indicate those 
with /cgraZ/pc asil = l-Oi 0-7, and 0.4. The constraints from the 
reddening E{B — V)/Nii are shown for /c,gra = 0.4, 0.7, and 
1.0 (with /pcasil = l-0)i drawn by thin (brown, green, and blue) 
contours. 

4.3. Inclusion of other carbonaceous grains 

Various populations of carbonaceous grains may con- 
stitute cosmic dust in addition to graphite. We explore 
the possible significance of other species of carbonaceous 
grains, such as glassy carbon or amorphous carbon. We 
find, however, that the two-component dust model, com- 
posed of glassy or amorphous carbon and astronomical 
silicate, does not give the extinction curve that lies in 
the 1 a ranges of observation for any /c.gic//Fo,asii and 
/c,amc//Fe,asii ratios. Thcsc carbouaccous grains do not 
give the proper UV bump at 2175 A. Glassy carbon shows 
a broad hump at around 2000 A, and amorphous car- 
bon has only a broad maximum at around 2000 A. We 
consider in the following the model in which graphite is 
partly replaced with glassy or amorphous carbon. 

Figure 7 shows a few example extinction curves for 
three-component models, compared with our fiducial 
graphite-silicate model, with q = 3.5 and amax ~ 0.24 
^m, fixed at the same parameters. While some of the 
resulting curves displayed do not satisfy 1 a range, as 
we replace graphite with these components excessively 
and/or we do not tweak the other parameters, we main- 
tain this set of parameters to see how the extinction curve 
is modified with the inclusion of the other carbonaceous 
material. 

Varying the (g, Omax) parameters, we examine to how 
much extent graphite can be replaced with glassy or 
amorphous carbon in the graphite-silicate model to main- 
tain the extinction curve within the 1 a range. Fig- 
ure 8 presents the maximum fraction of C atoms in 
glassy (/cgic) or amorphous carbon (/c.amc) relative 
to the entire condensed component of carbon /c. grain 

against /c,grain/./Fc,asil, where /c.grain = /cglC + /Cgra 
or /Cgrain = /c,amC + /c,gra- Wc SCC that Up tO «30- 

40 % or R:!50-60 % of carbon in grains can be in glassy 
or amorphous phase, respectively. In other words, more 
than 60-70 % or 40-50 % of C atoms must be in graphite. 
With this inclusion, the parameters (g, flmax) are shifted 
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Fig. 7. — (a) Extinction curves from the graphite-glassy (or 
amorphous) carbon-siUcate models with q = 3.5 and amax = 
0.24 fim, and (b) the expanded figure of (a) for NIR, as com- 
pared with our fiducial graphite-silicate model (thick black curve). 
/c,grain//Fc,asil = 0.7 is taken for all cases. The 1 cr ranges from 
the observed extinction curves in FM07 are indicated with error 
bars. 

only a little; see Figure 9. Our representative parame- 
ters q = 3.5 and Omax = 0.24 fim still remain to be in the 
allowed solution. 

Silicon carbide (SiC) is another candidate material 
that can be in dust grains. We find the upper limit 
of carbon contained in SiC being approximately 15 % 
to reproduce the extinction curve (see Figure 8). The 
abundance of SiC is also limited from the availability of 
Si, which is 11 % of the abundance of C. We note that 
the constraint from the extinction is similar to the limit 
from the abundance ratio. Graphite must be the major 
condensation of carbon. From the point of view of Si, on 
the other hand, a large fraction of Si can be in SiC with- 
out disturbing much the extinction curve: /si,sic can be 
unity for /^grains = 0.61-0.87. The (g, Cmax) parame- 
ter range that gives 1 cr of extinction becomes slightly 
shrunk but basically unchanged (see Figure 9). 

With SiC included, however, the parameter region al- 
lowed from E{B — V)/Nii is appreciably squeezed, leav- 
ing the overlap from the size of E{B — V)/N-ii and the 
extinction curve only at the lower edge of the (q, flmax)- 
g = 3.5 and Omax = 0.24 fim are not a solution to sat- 
isfy the two constraints any more. The inclusion of SiC 
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Fig. 8. — Allowed fractions of C atoms in the glassy phase 
(/c,gic//c, grain, blue), amorphous phase (/c,amc//c, grain, red), 
and silicon carbide (/c,Sic//c, grain, green). The abscissa is 

/C,grain//Fe,asili where /c, grain = /c,gIC + /c,gra Or /c, grain = 

/c ,amC + /Cgra or fc , grain /cSiC + /c,gra- The ordinate is the 
maximally allowed fraction of the individual component. 
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Fig. 9. — Allowed regions of q and Omax with which the 
model satisfies the 1 a extinction ranges for the MW when glassy 
carbon, amorphous carbon, or silicon carbide is included in the 
graphite-silicate model. The parameters are marginalised over 
/c,grain//FG,asil- The regions that satisfy the observed size of red- 
dening E{B ~ V) are also shown (dashed curves) for /c, grain = 0.7 
and /pc asil = l-O. The blue, red, and green lines are, respectively, 
for the cases where glassy carbon, amorphous carbon, and SiC are 
added to the graphite-silicate model, which is represented with the 
black curves for comparison. For the E(B — V)/Nji constraint it 
is taken that /c.gic//c,grain = 0.3, /c,amc//c,grain = 0.3, and 

/c,Sic/.fc, grain = 0.1. 



necessarily reduces the graphite and silicate components, 
whereas the Q factor of SiC is significantly smaller, and 
SiC contributes little to the cross section. 

In conclusion, the inclusion of amorphous or glassy car- 
bon disturb little the agreement of the graphite-silicate 
model with the observation. The inclusion, however, 
does not make the overlap of the two constraints easier. 
The dust to gas ratio is unchanged from the graphite- 
silicate model. With SiC the overlap of the two con- 
straints becomes marginal, and the dust to gas ratio be- 



comes Pdust/PH = 1/(122+^°). 
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Fig. 10. — Central wavelength and the width parameter that 
describe the 2175 A feature of the MW extinction. The observed 
data from FM07 (squares) are compared to our model with the 
parameters that gives a good fit (within 1 ct of the full extinction 
curve) to the global extinction curve. The horizontal and vertical 
lines indicate 1 ct of the data. 

4.4. The 2175 A Feature 

With small graphite the feature at 2175 A is generated. 
It is known that this feature is observationally fit well 
with the Drude formula, including a smooth background, 



Av 



1 

Ry 



co + ciA + 



C2 



(13) 

In Figure 10 we give the central wavelength xo = Aq ^ 
and the width of the profile 71 for 328 extinction curves 
of FM07, which are summarised as 



Ao = 2178±10A 71 = 0.95 ±0.1 /im" 



(14) 



We fit the feature from our graphite-silicate model, which 
reads 



An = 2186 ± 12 A 



7i = 1.21 ±0.08 fim- 



(15) 



indicating that the model width is broader by 30 % (see 
also Draine & Malhotra 1993), while the central wave- 
length agrees with the observation. We also show the 
model with amorphous or glassy carbon included with its 
fraction and size parameters in the range allowed from 
the entire extinction curve. Their inclusions make the 
width slightly (wlO %) smaller, but not sufficient to give 
the observed width. This seems to be a problem intrin- 
sic to the optical data we adopted for graphite. More 
detailed treatment may be needed for small graphite or 
PAH. 

4.5. Inclusion of Fe, Fe^Oi, FeS, and AI2O3 

We consider how much iron can be incorporated in 
generic iron-bearing grains rather than in astronomi- 
cal silicate. Figure 11 shows the extinction curve for 
the model with our typical parameters q = 3.5 and 
flmax = 0.24 fim. Figure 12 exhibits the upper limits 
of the fraction of Fe in metallic phase (/Fe,Fe, solid line) 
and in magnetite, Fe304 (/fo,Fo304i dashed line), for the 
model that satisfies the 1 cr ranges of the extinction. 
The rest of Fe atoms in interstellar space are assumed 
to be condensed in astronomical silicate, i.e., /fc, grain = 

/Fe.Fe + /Fc,asil = 1 Or /fo. grain = /fc,Fo304 + /Fc.asil = 1, 
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Fig. 11. — (a) Extinction curves for the graphite-silicate model 
where iron-bearing silicate is replaced with metallic Fe, Fe304 or 
FeS, while the corresponding silicate portions are replaced with 
forsterite (Mg2Si04). We take the grain-size parameters, q = 3.5 
and amax = 0.24 /^m. (b) the expanded figure for the NIR. We 

take /Cgrain/Zpo, grain = 0.7. 

and Mg not included in astronomical silicate is assumed 



to be in forsterite, so that ^ 



Mg 



1. Here, the 



abscissa is chosen to be the C/Fe ratio of the condensa- 
tion. 

There are parameters (/c,gra//Fc, grain = 0.15-1.43), for 
which Fe may stay entirely in the metallic phase rather 
than in astronomical silicate, yet the model gives the 
correct extinction curve. A similar parameter range, 
however, is very small with the Fe304-graphite- forsterite 
model. Practically, the maximum abundance of FC3O4 
is about 80 %. The NIR extinction of Fe3 04 is large 
by more than a factor of 2 compared with the observed 
extinction, whereas Fe grains give the extinction in NIR 
only by «30 % larger, and, therefore, the abundance of 
magnetite is limited from the extinction curve in NIR 
more strongly. The 1 cr-allowed (g, a,nax) parameters 
with these models are shown in Figure 13. 

The combination of metallic Fe and forsterite works in 
a way virtually the same as astronomical silicate, which 
means that the Mg/Fe ratio in astronomical silicate is ar- 
bitrary. The situation is somewhat different if iron is in 
magnetite. The inclusion of Fe304 disturbs the allowed 
range in ((7, flmax) plane by an appreciable amount, mak- 
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TABLE 3 
Allowed Grain Species 



Dust Grains 


1 M07 1 cr 


CCM 1 a 


(1) Graphite— Astronomical Silicate 


Yes 


No 


(2) Glassy Carbon — Astronomical Silicate 


No 


No 


(3) Amorphous Carbon— Astronomical Silicate 


No 


No 


(4) Graphite-Glassy Carbon-Astronomical Silicate 


Yes 


No 


(5) Graphite-Amorphous Carbon-Astronomical Silicate 


Yes 


No 


(6) Graphite-SiC-Astronomical Silicate 


Yes 


No 


(7) Graphite-Fe 


Yes 


No 


(8) Graphite- Fes O4 


No 


No 


(9) Graphitc-Fc— Astronomical Silicate 


Yes 


No 


(10) Graphitc-Fc304-Astronomical Silicate 


Yes 


Yes 


(11) Graphit0-Fe-Mg2SiO4 


Yes 


No 


(12) Graphite-Fes 04-Mg2Si04 


Yes 


Yes 


(13) Graphite-Fe— Astronomical Silicate-Mg2Si04 


Yes 


No 


(14) Graphite-Fes 04-Astronomical Silicate-Mg2Si04 


Yes 


Yes 


(15) Graphite-FeS-Astronomical Silicate-Mg2Si04 


Yes 


No 



9) 



Note. — We count the species when its fraction is more than 10 % 

0.4 
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f /f 

C.gra' Fe, grain 



Fig. 12. — Maximum fractions of Fe atoms allowed in the 
metallic (/pe.Fe //pe, grain i red solid) or the magnetite phase, Fe304 
(/PG,PG3 04//pe,graini blue dashed), with which the model satis- 
fies the 1 (T extinction range for the MW. The maximum frac- 
tion is also drawn for FeS (/po,Pos//po,graini green dotted). The 
abscissa is /c,gra//Pe jgrain ; 

where /pe 

jgrain — /Pe.Pe + /pe , asil 5 
/Po, grain = /Po.Pog O4 + /Pe,asil i Or /Pe, grain = /Pe.PeS + /Pe.asih and 

the ordinate is the fraction of each species added to the model. 

ing flmax smaller by 20 %. Nevertheless, we see that the 
overlap of the two constraints is maintained. 

The cosmic abundance of sulfur is about 0.7 times that 
of Fe. We consider the case where all S is in troilitc, FcS. 
The rest of Fe is in astronomical silicate and further the 
rest of Mg is in forsterite. Troilite also gives a large NIR 
extinction almost as much as magnetite, but maximal 
amount of FeS is allowed because of smaller cosmic abun- 
dance of sulfur. This case satisfies the 1 cr constraint with 
flmax somewhat smaller than the graphite-silicate model, 
as seen in Figure 13. 

In conclusion iron can be in a variety of grain species, 
including olivine, metallic phase, Fe304, or FeS without 
disturbing the extinction curve. The only condition is 
that Fe is not predominantly in magnetite, which pro- 
duces too large NIR extinction. The grain size parame- 
ters are nearly the same for all cases, up to the result that 
the maximum size cutoff becomes by 20-30 % smaller if 
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Fig. 13. — Allowed regions of q and amax with which the 
model satisfies the 1 cr extinction ranges for the MW, marginalised 

over /c,gra//Pc, grain- We take /c.gra = 0.7 and /pc, grain = 1-0. 

The blue, red, and green curves are, respectively, for the inclu- 
sion of metallic Fe, magnetite (/PC3O4//PC, grain = 0-8), and FeS 
(/pcs//pc, grain = 0.68). The black curve is our fiducial graphite- 
silicate model. 

Fe304 or FeS is the major component. All iron atoms 
need not necessarily be locked in astronomical silicate, 
and the ratio of Fe:Mg in astronomical silicate is ar- 
bitrary. Allowing for the inclusion of a variety of iron 
material, the mass density of dust differs little from the 
graphite-silicate model. 

We can ignore the contribution from corundum 
(AI2O3). In addition to the small abundance of Al (1/13 
of Mg), the Q factor for corundum is small (Toon et al. 
1976). Even if all Al atoms are locked in corundum, they 
contribute little to modifying the extinction curve. 

Table 3 summarises our result for most grain species 
wc considered in the present work. Limiting ourselves to 
at most four grain species, we consider what combina- 
tion of species would give the extinction curve consistent 
with the observation at 1 cr. Some of them were already 
discussed in the text above. We take the species as valid 
when its fraction is more than 10 %. We here note one 
particular case that can reproduce the extinction law of 
CCM: it is the combination of graphite, astronomical sili- 
cate, and magnetite. This is due to the particularly large 
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Fig. 14. — (a) Extinction curves for the SMC from the graphite- 
silicate dust models with q = 3.5 and amax = 0.24 /xm, and (b) 
the expanded figure of (a) for the NIR. The solid, dashed, and 
dotted curves are, respectively, for /c,gra//Fc asil = Oj 0.1, and 0.4 
with ajnin = 0.005 [ira fixed as our fiducial choice. The dot-dashed 
line is with the minimum cutoff increased to amin = 0.02 /xm, 
while /c,gra//Fe,asil = 0.4. The error bars span the maximum and 
minimum (taken as '1 cr') of the observed SMC extinction curves 
of G03. 

Q factor of magnetite in the NIR, whieh is excluded if 
we take FM07. 

5. DUST IN THE SMC 

Figure 14 shows examples of the model extinction 
curve for the SMC with q = 3.5 and Omax = 0.24 
/xm for /c,gra//Fc,asii = 0, 0.1, and 0.4. (The curve 
for /c,gra//Fo,asii = 0.4 appears to be out of the 1 a 
range, but satisfies it if Omax — 0.19 [ivn.) The model 
for /c,gra//Fo.asii < 0.4 rcproduccs the SMC extinction 
curve within 1 a when q and Omax are adjusted. With 
/c,gra//Fe,asii = 0.4 One sces in this figure the symp- 
tom that a small 2175 A bump starts appearing and the 
rise towards the far UV side becomes insufficient. We 
draw another curve (dot-dashed curve) with the mini- 
mum grain size cutoff increased to Omin = 0.02 /im, re- 
moving much of small graphite grains. The bump at 2175 
A is reduced and becomes insignificant, but the removal 
of small size grains makes the far UV rise insufficient. 
Therefore, such a case is excluded. Within the graphite- 
silicate model, the way to make the predicted extinc- 
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Fig. 15. — The abundance of C and Si in dust grains relative 
to hydrogen for the SMC to meet the extinction curve. The range 
left to the oblique line is allowed from the extinction curve, and 
that within the two curves running nearly horizontal is from E(B — 
V)/Nii. The two dashed lines (horizontal and vertical) show the 
total abundance of C and Si from Russel & Dopita (1992). 

tion consistent with the observation is not to change the 
grain size parameters but to decrease the abundance of 
graphite grains. The steeper rise towards far UV is ac- 
counted for by the dominance of silicate grains in their 
small size. 

The abundance of Si and C contained in dust is shown 
in Figure 15, where the abundance from Russell & Do- 
pita (1992) is also indicated with dashed lines. We have 
log(Si/H) + 12 = 6.7-7.8 and log(C/H) + 12 < 7.9 
from the extinction in the SMC. When /si, grain = 1, 
log(C/II) + 12 < 7.6 in agreement with eq. ([TS]) be- 
low. It is noted that the carbon abundance in graphite 
in the SMC is at least by a factor 1.5 times less than 
the abundance of Russell & Dopita (if Si is all condensed 
into dust; following the horizontal dashed line); the lat- 
ter gives 7.73 for C and 7.03 for Si. The upper and lower 
curves correspond to the 1 a error of N-h_/E{B — V). For 

the SMC Pgra/Pasil = 0.54/c.gra//Fo,asil < 0.22. 

Figure 16 is a summary of the region of q and a^aax 
with which the graphite-silicate model gives the extinc- 
tion curve within the 1 cr extinction range marginalised 
over /c,gra//Fo,asii- Thc allowed grain size lies again for 
3-3 £ g < 3.8 and 0.19 /.tm < Cmax ^ 0.35 /im. It is 
interesting to observe that this range for the SMC agrees 
with that for the MW in spite of the significantly differ- 
ent behaviour of the two extinction curves, as was noted 
earlier in Pei (1992). 

The other curves in Figure 16 show the regions allowed 
by Nii/E{B - V) for the 1 a error of Eq (2). We find an 
overlap of the two constraints for 



/Cgra < 0.41. 



(16) 



This maximum allowed fraction of graphite is smaller 
than thc minimum /c,gra//Fe,asii favoured for the MW 
(/cgra > 0.56), suggesting that condensation of carbon 
into grain should be less efficient in the SMC, as might 
be satisfied if /c,gra oc Nq- More detailed arguments, 
however, depend on the accuracy of the abundance esti- 
mates, in particular, for the SMC. 

For the model that satisfies 1 a constraints of both 
extinction curve and E{B — V) size for the SMC, we find 
Pdust/pH = 1/760 ]'^ggQ, where the large error range comes 
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Fig. 16. — Allowed regions of q and amax with which the model 
satisfies the '1 cr' extinction ranges for the SMC. The graphite- 
to-silicatc ratios is marginalised over < /c,gra//Fe,asil ^ 0.41. 
The inlaid curves indicate the regions for /c,gra = 0, 0.2, and 
0.4. The regions for 0.4 is a point indicated with a small circle 
(the actual region is enlarged) near the bottom of the allowed re- 
gion. The constraints from the reddening E{B — V)/Nii are shown 
with brown, green, and blue contours for /c,gra = 0, 0.2, and 0.4 
(/pc asil = 1.0). The thin dotted curve recapitulates the allowed 
region for MW dust in Figure 6 for comparison. 

from that of E{B — V)/N}i, and corresponds roughly 
to the upper and lower curves of the allowed region in 
Figure 15. If the abundance is constrained with Russell 
& Dopita's value and Fe is all condensed, 

+70 

Pdust/PH = 1/760 (17) 

We have 

iCext = (2.2 ± 0.3) X lO'^ mag cm^ g'^ (18) 

(It is (2.1 ± 0.4) X lO"' mag cm^ g~^ if the abundance 
constraints are removed.) The ratio of graphite to sili- 
cate is Pgra/Pasil = 0. 1 1 ± 0. 1 1 , Or Nc/Npc = 1.6 ±1.6, in 
numbers of atoms. This ratio is 4.5 times smaller than 
that for the MW. The abundance itself indicates the ra- 
tio of C/Fe to be 1.4 times smaller, so the formation of 
carbonaceous grain is suppressed by ?a3 times more in 
the SMC. 

6. CONCLUSIONS AND DISCUSSION 

We have confirmed that the graphite-silicate grain 
model gives a satisfactory description of the extinction 
curves within the simple power-law model of the size 
distribution. The grain size distribution is tightly con- 
strained to the index q = 3.5 ± 0.2. We showed that 
departures from the power law is not needed. However, 
we necessarily need the cutoff of the power-law at some 
maximum size a = 0.2-0.3 fiui. Grains may be of a va- 
riety of species of carbonaceous and silicate materials. 
Their size parameters, however, vary little and the same 
parameters also apply to both MW and SMC extinction 
curves in spite of their apparently different behaviours. 
The difference of the two extinction curves seems to be 
ascribed to the abundance of graphite relative to silicates, 
and hence to some lower efficiency (< 1/2) of graphite 
condensation further to the lower ratio (< 1/1-5) of car- 
bon to silicate in the SMC indicated in the abundance 



& Fukugita 

estimate currently available. We remarked the somewhat 
dissatisfying description of the 2175 A feature. 

While grains can be a variety of combinations that 
contain silicate and carbonaceous material, the presence 
of a significant abundance of graphite seems important 
in any case. It is interesting to note that the resulting 
q ~ 3.5 is the power expected from coUisional equilibrium 
for small grains. We have derived elemental abundance 
of Si and C contained in grains from the extinction. It 
takes for the MW values consistent with the solar. If 
we take the widely adopted abundance of GS98, we infer 
the fraction of carbonaceous grains against total carbon, 
/Cgra = 0.6-0.7, which agrees with the observational de- 
pression factor for carbon. This may be compared with 
the corresponding value for the SMC /c,gra S 0-4 tak- 
ing the elemental abundance of Russell & Dopita (1992). 
Required silicate is also consistent. 

The upper cutoff is compelling for the grain size distri- 
bution to give the correct shape of the extinction curve. 
Extending it to a larger size would largely dismiss the 
agreement with the observation (sec Drainc 2009): for 
example, the Ry parameter becomes intolerably large. 
While we are not able to find the reason for the cutoff, we 
may see another argument that forces us to impose the 
presence of a cutoff in the power law. The abundance of 
grains is constrained at a submicron region to give cor- 
rectly the size of the observed optical extinction. The 
power law with q = 3.5 means that the mass of grains is 
slowly increasing with the upper cutoff mass rric of the 
integral as m°'^^. If the cutoff were larger, we would have 
too large a mass in grains to be accounted for as a prod- 
uct of stellar evolution. The dust mass density relative 
to hydrogen Pdust/pn ~ 0.008 means fidust ~ 4.4 x 10~^ 
taking the global hydrogen abundance for HI and H2 ob- 
servations, fin = 5.4 X 10~*. If dust is a product of stel- 
lar evolution over the cosmic time, stars of the amount 
r^star — 0.003 would producc dust grains no more than 
i^dust ~ 1 X 10~^ (Fukugita 2011). This amount is con- 
sistent with the power-law distribution constrained from 
the extinction if the power-law is cutoff at 0.25 /xm. If 
the maximum size cutoff were 3 times larger, say, the 
dust abundance would be larger by ~2 times, more than 
the star formation activity can account for. 

Small astronomical objects may also obey the same 
power law. Their normalisation, however, should be 
smaller by a large factor than that for the dust grains. 
The integrated mass density of the core of planets for 
a > 100 km is estimated to be 1/300 that of dust grains 
(Fukugita 2011). The power law should break by two 
orders of magnitude at a few tenth of micron. The addi- 
tion of this planet core mass disturbs little the estimate 
of the mass density of grains. 

We have discussed that the behaviour of extinction 
curve in NIR is important. If the power of NIR extinc- 
tion curve A"'' is as small as 7 = 1.6, as was derived 
earlier and taken by CCM or WDOl, one cannot repro- 
duce the extinction curve from UV to NIR by the grain 
model with the simple power law of the grain size. One 
need to devise a substantial tweaking of the size distri- 
bution from, or extra components added to, the power 

This is caused by the fact that the Ry parameter resulting 
from silicate grains is wildly oscillating as a function of a for a > 0.2 
fim, which makes the resulting Ry parameter too large. 
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law as was done by WDOl, Zubko et al. (2004), and oth- 
ers. With a larger power of the NIR wavelength depen- 
dence, such as 7 = 1.8-2.0, however, the simple power 
law model works to account for the observed extinction. 
The NIR power 7 = 1.6 could be consistent with the 
power law grain size, only when the iron component of 
grains would largely be in magnetite, which has a large 
NIR scattering efficiency. 

We have seen that the presence of graphite grains is 
uniquely important, at least in the MW, although the 
condensation of C atoms may not necessarily be all into 
graphite. Roughly half amount of carbon may be con- 
densed into a glassy or an amorphous phase, or by a 
non-negligible fraction in silicon carbide, whereas some 
significant fraction of C atoms must be in graphite. Iron 
may also remain in the metallic phase. There is no need 
for all Fe atoms to be locked in astronomical silicate. 
The effective ratio Mg to Fe in olivine seems arbitrary. 



Fe304, however, cannot be predominant. Fe may also be 
in troilite, as much as the sulfur abundance allows. 

The size distribution is well converged to a narrow 
range, whatever are other grain compositions included. 
The variation of extinction along lines of sight may be 
accounted for by a small variation of the grain-size dis- 
tribution with Aq « 0.2 or Aamax/amax ~ 0.3 and with 
the graphite to silicate ratio as is the difference between 
the MW and the SMC. It is noteworthy that the a~^-^ 
power as expected in coUisional equilibrium seems to be 
generic to dust grains. 
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